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ABSTRACT 

The Hubble Advanced Camera for Surveys has been used to obtain a narrow-band image of the 
weak emission peak seen at A = 7205 A in the Gunn-Peterson Ly/3 absorption trough of the highest 
redshift quasar, SDSS J1148+5251. The emission looks perfectly point-like; there is no evidence for 
the intervening galaxy that we previously suggested might be contaminating the quasar spectrum. We 
derive a more accurate astrometric position for the quasar in the two filters and see no indication of 
gravitational lensing. We conclude that the light in the Ly/3 trough is leaking through two unusually 
transparent, overlapping windows in the IGM absorption, one in the Ly/3 forest at z ~ 6 and one in the 
Lya forest at z ~ 5. 

If there are significant optical depth variations on velocity scales small compared with our spectral 
resolution (~ 150 km/s), the Lya trough becomes more transparent for a given Ly/3 optical depth. 
qq ■ Such variations can only strengthen our conclusion that the fraction of neutral hydrogen in the IGM 

increases dramatically at z > 6. We argue that the transmission in the Ly/3 trough is not only a more 
sensitive measure of the neutral fraction than is Lya, it also provides a less biased estimator of the 
' neutral hydrogen fraction than does the Ly a transmission. 

Subject headings: comology: observations — early universe — galaxies: high-redshift — intergalactic 
medium — quasars: absorption lines — quasars: individual (SDSS J1148+5251) 

1. INTRODUCTION 

The epoch of reionization of intergalactic hydrogen is a key cosmological observable that constrains the star and galaxy 
' formation history of the universe. Clues to the ionization state of the intergalactic medium (IGM) come from two principal 
^ , sources: cosmic microwave background (CMB) polarization measurements, which are sensitive to the integrated electron 
O ,■ scattering optical depth from z = to z ~ 1000, and studies of hydrogen absorption in the spectra of quasars and galaxies, 
L_ ' which can detect neutral gas along the line of sight at z < 6. 

The situation is currently somewhat muddled. The CMB measurements from the first year of Wilkinson Microwave 
Anisotropy Probe (WMAP) observations indicate that the IGM may already have been ionized at z > 10 (Kogut et al. 
^ | 2003). Quasar spectra show evidence for a sharp increase in the Lya absorption at z > 6, leading to the suggestion that 
reionization occurred near that time (Becker et al. 2001; Fan et al. 2002, 2003, 2004; White, Becker, Fan & Strauss 2003; 
but see Songaila 2004 for an opposing view). Wyithe & Loeb (2004) and Mesinger & Haiman (2004) conclude from the 
extragalactic H II regions around these quasars that the IGM is indeed close to neutral. On the other hand, the mere 
jJJ ■ detection of Lya-emitting galaxies at higher redshifts is argued to be inconsistent with the IGM having a significant 
d \ neutral fraction at z ~ 6.5 (Rhoads et al. 2004; Stern et al. 2004; Malhotra & Rhoads 2004). 

Theoretical work has shown that the WMAP and quasar results are not necessarily contradictory because they are 
sensitive to rather different physical regimes (Gnedin 2004). The black Gunn-Peterson (1965) Lya absorption troughs in 
quasar spectra can be produced by a small fraction (~ 10 -2 ) of neutral hydrogen, while the CMB polarization can be 
attributed either to an interval at high redshift when the universe was ionized followed by another recombination, or to a 
more recent reionization that persisted to the present day (see Fig. 7 of Gnedin 2004.) There are also theoretical reasons 
to believe that the ionization history of the universe could have been complex, with multiple periods of ionization and 
recombination (Wyithe & Loeb 2003, Cen 2003). Recent studies indicate that even a completely neutral IGM would not 
necessarily suppress Lya emission from galaxies at z ~ 7 (Haiman 2002, Santos 2004, Gnedin & Prada 2004), especially 
when such uncertain effects as galactic winds, clustering of star-forming galaxies, and radiative transfer in the scattered 
Ly a radiation are included. We need better observations to distinguish among the many possibilities. 

Our understanding of quasar absorption at z > 6 is based on only a handful of quasars, so cosmic variance along the 
different lines of sight is clearly a concern. There are just four quasars known with redshifts beyond 6.2, all discovered 
using the Sloan Digital Sky Survey (Fan et al. 2004), and one of those has intrinsic broad absorption lines that make 
it unsuitable for studies of the intervening Lya absorption. White et al. (2003) presented high resolution, high signal- 

1 Based on observations obtained with the NASA/ESA Hubble Space Telescope (HST). HST is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under NASA contract NAS5-26555. 

2 Space Telescope Science Institute, Baltimore, MD 21218 

3 Physics Dept., University of California, Davis, CA 95616 

4 IGPP/Lawrence Livermore National Laboratory 

5 Steward Observatory, University of Arizona, Tucson, AZ 85721 

6 Princeton University Observatory, Princeton, NJ 08544 
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to-noise ratio Keck spectra of the two highest redshift quasars, SDSSp J103027. 10+052455.0 (z — 6.28, hereafter SDSS 
J1030+0524) and SDSS J114816.64+525150.3 (z = 6.42, hereafter SDSS J1148+5251). SDSS J1030+0524 shows black 
hydrogen absorption with no indication of any residual quasar light in either the Ly a or Ly (3 troughs over the redshift 
range 5.97 < z < 6.18, which sets strong limits on r Q , the Lya optical depth derived from the mean transmission. This 
spectrum has been analyzed by several different groups (e.g., Mesinger & Haiman 2004) and is the strongest evidence that 
the IGM is becoming neutral at z > 6. 

The spectrum of SDSS J1148+5251 also shows deep Lya and Ly (3 absorption but differs from SDSS J1030+0524 in 
that the absorption troughs are not completely black (Fig. 1). There is a weak peak at z = 6.08 that is seen in both 
Ly a and Ly (3 and is likely due to an ionized bubble along the line-of-sight that produces a "leak" in the IGM absorption 
(White et al. 2003). There are also several stronger emission peaks in the Ly/3 trough without associated Lya peaks, 
with the strongest at 7205 A. The straightforward interpretation of these peaks as due to transparent Ly (3 windows in 
the IGM faces some serious problems, however (discussed further in §4.) Consequently, White et al. suggested that these 
features might be Ly a emission at z ~ 5 from an intervening galaxy; the presence of strong C IV absorption at the same 
redshift provided some support for this suggestion. The implausibility of finding such a galaxy along the line-of-sight by 
chance would be reduced if the intervening object is lensing the quasar. The uncertainty about the nature of the residual 
light in the Lya and Ly/3 troughs of SDSS J1148+5251 called into question its use for studying the ionization state of 
the IGM at z ~ 6. 

In this paper, we report Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) observations of the quasar 
SDSS Jl 148+5251 using a narrow filter in the Ly/3 trough centered on the 7205 A peak. These ACS images discriminate 
between the two hypotheses by allowing us to look directly for evidence that the emission comes from an intervening 
galaxy and to check for any signs that the quasar is gravitationally lensed. The observations are described in §2, and §3 
discusses the analysis. The results and their implications for this object and for the IGM ionization are dicussed in §4. 

2. OBSERVATIONS 

Images of SDSS J1148+5251 were obtained with the ACS Wide Field Camera (WFC) on 2003 October 25. Two ramp 
filters were used: a narrow band covering the emission peak at 7205 A (FR716N-7220, 4256 s exposure time) and a broader 
band located on the quasar's Lya emission line (FR914M-9050, 720 s). We briefly refer to these filters as F720N and 
F905M in the following. Figure 1 shows the location of the filter bandpasses with respect to the quasar's spectral features. 
About 60% of the F720N counts come from the narrow feature at 7205 A, with most of the remainder coming from the 
weaker features nearby. The short-exposure F905M image is used to determine accurately the position and morphology 
of the quasar for comparison with the F720N detection. 

Note that the field of view in the ACS ramp filters is limited to about 40 x 80 arcsec. The filter strip is narrow, covering 
only about a third of the ACS field of view, and the ramp central wavelength changes along the length of the filter so 
that different parts of the image are being observed at different wavelengths. This does complicate the data analysis, but 
the limited field does not interfere with our science goals. 

The images have been processed using the ACS Science Team's APSIS pipeline (Blakeslee et al. 2003), which aligns 
the images, corrects the astrometry to the GSC-2 coordinate system, produces combined drizzled images on the same 
output pixel grid for the two filters, and runs SExtractor (Bertin & Arnouts 1996) to identify sources in a detection image 
created by combining the two filter images. We have augmented the processing using a custom cosmic-ray (CR) rejection 
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Fig. 1.— Keck ESI spectrum of SDSS J1148+5251 (White et al. 2003) with the filter throughputs for our HST ACS observations. The narrow 
blue filter (FR716N-7220) is placed in the dark Ly/3 absorption region around z = 6 and isolates the narrow emission peak at A = 7205 A. 
The broader red filter (FR914M-9050) encompasses the quasar's strong Lya emission line and continuum and is used as a reference image. 
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algorithm that is optimized for the case of only two CR rejection images. The algorithm recognizes instances where cosmic 
rays hit the same pixel in both images and flags both as bad. The observations included two CR-split exposures at each 
of two different dither positions for F720N and a single exposure at the two dither positions for F905M. 

Figure 2 shows the resulting ACS images. The large image is the detection image while the insets show the quasar in 
the two filters. The detection image is dominated by the long-exposure F720N image except for the quasar, which is much 
brighter in the F905M image. The count rate for the quasar is 0.83 counts/s in F720N and 54.5 counts/s in F905M. 

An inspection reveals no indication that either the F720N or F905M image is resolved, nor is there a significant shift 
between them. We quantify this statement in the next section, but our basic conclusion is that the 7205 A emission peak 
is perfectly consistent with quasar light leaking through a transparent window in the z ~ 6 Ly/3 and z ~ 5 Lya IGM 
absorption. 

3. ANALYSIS 

Our principal goal in acquiring these observations was to determine the nature of the emission seen in the Ly (3 absorption 
trough at 7200 A. We therefore focus on determining whether the emission is extended (as expected for an intervening 
galaxy) or point-like (as expected for quasar light) and whether there is any indication of an astrometric shift between 
the two images (which would be expected if the quasar were gravitationally lensed.) We also determine accurate absolute 
astrometry for the quasar using nearby SDSS objects as references. 

3.1. Size Limits 

An upper limit on the size of the F720N image was determined by fitting a monochromatic 7205 A Tiny Tim (Krist 
1995) point-spread function (PSF) to the F720N quasar image. The simulated PSF was subsamplcd 7 x 7 to allow an 
accurate sub-pixel position match. The fitting was performed on the flat-fielded image rather than the drizzled image to 
avoid possible problems with the drizzled PSF and to incorporate the best possible knowledge of the noise characteristics 
of the data. (Fits to the drizzled image gave similar results.) 

The data are fitted very well by the PSF plus sky; the reduced y 2 for a 20 x 20 pixel region (lxl arcscc 2 ) around the 
quasar is less than unity. To determine the size limit, the PSF was blurred with a Gaussian. When the Gaussian width 
is too large, it is no longer possible to get a good fit as determined by the x 2 value. The 90% confidence upper limit on 
the source size is a FWHM of 0.059 arcsec, and the 99% confidence limit is 0.067 arcsec. 

For comparison, galaxies at z — 5 typically have half-light diameters of 0.5 arcsec, and are almost never seen to have 




Fig. 2. — HST ACS/WFC images of SDSS J1148+5251. The large image shows the detection image (a combination of the two filters), 
while the insets depict the quasar image in the two filters. The quasar appears unresolved, and the images coincide to high accuracy. 



4 



White, Becker, Fan, Strauss 



diameters smaller than 0.2 arcsec (Ferguson et al. 2004). Our 99% confidence size limit is far smaller than would be 
expected for an intervening z — 5 galaxy. 

While we can confidently conclude that most of the light in F720N does not come from a galaxy, we cannot rule out 
the possibility that a small fraction of the light is extended. Fitting the sum of a point source and a 0.5 arcsec FWHM 
Gaussian reveals that even a galaxy with 30% of the total flux would not be detectable. The 90% confidence upper limit on 

the fraction of light contributed by a galaxy is 0.46. That corresponds to a continuum level of 2.2 x 10~ 19 erg cm~ 2 s -1 A . 
The best limit on residual continuum emission from an intervening galaxy probably comes from the spectrum itself and 
is about 5 x 10~ 20 . 

3.2. Relative Astrometry 

The relative astrometric alignment of the two filter images was established by comparing positions for all common 
sources in the field of view. We fitted Tiny Tim PSFs to the drizzled images to determine the quasar position in each 
filter. (In this case we used the drizzled images because we needed to correct for the ACS distortion.) The measured 
quasar positions differ by only 0.22 pixels (0.011 arcsec). This is very similar to the rms for the objects matched to align 
the images. We conclude that the images in the two filters have positions that agree to better than 0.02 arcsec. 

Note that if a lensing galaxy produced some of the light being detected in the F720N filter, we would expect to see some 
apparent shift between the positions because the apparent lensed quasar position would be shifted away from the galaxy 
center, and the galaxy would be far more prominent in the F720N filter. There is no hint of any such effect. Richards et 
al. (2004) also found no evidence for lensing in HST ACS images of four other high-redshift quasars. 

3.3. Absolute Astrometry 

We corrected the absolute astrometry of the HST images by matching our SExtractor catalog with the SDSS DR3 
photometric catalog (Abazajian et al. 2004) and measuring the mean shift between them. We excluded the quasar itself 
from this matched list to avoid biases. There are 29 ACS/SDSS matches within 50 arcsec of the quasar. We searched for 
matches with separations as large as 2 arcsec, but all the SDSS sources found fell within 0.6 arcsec of the (APSIS-corrected) 
ACS positions. We conclude that all the matches are reliable associations. 

The coordinate offsets in RA and Dec were estimated using both mean and median statistics, with essentially identical 
results. Both RA and Dec offsets are accurate to 0.03 arcsec, which is our remaining systematic error compared with 
the local SDSS coordinate system. Of course, this does not mean that our absolute astrometry is that good, but the 
astrometric accuracy for the quasar position is about 0.1 arcsec, the limit set by the SDSS reference frame (Pier et al. 
2003). 

The corrected position for the quasar SDSS J1148+5251 is a = 11 48 16.645, S = +52 51 50.21 (J2000) with uncertainties 
of 0.003 s in RA and 0.03 arcsec in Dec. This differs by 0.13 arcsec from the quasar position given in the DR3 source 
catalog. Note that the extreme colors of this object (it is detected only in the z-band longward of 9000 A) introduce 
significant differential refraction in the ground-based SDSS images, which reduces the accuracy of the SDSS astrometry. 
Differential refraction is of course not an issue for HST observations. 

This improved position is most valuable in comparing to observations at other wavelengths. Walter et al. (2004) resolved 
the CO emission associated with SDSS J1148+5251 into two components. The CO positions are accurate, but the SDSS 
position falls between the components with uncertainties that overlap both. Walter et al. speculated that the quasar 
should be aligned with their southern component based on the overall symmetry of the source. Our improved position 
confirms that, with the quasar position falling exactly on top of the southern CO component. 

4. DISCUSSION 

Our conclusion from the ACS images is that there is no evidence for a galaxy along the line-of-sight to SDSS J1148+5251. 
The emission peak detected in the Gunn- Peterson Ly (3 trough at 7205 A is light from the quasar that is leaking through 
transparent regions in the IGM. 

We recognize that our arguments for an intervening galaxy were met with some (apparently well deserved) skepticism 
among our readers 7 . But the superficially more plausible explanation that the light is coming from the quasar actually 
faces some serious difficulties that were discussed in White et al. (2003). 

The principal difficulty is the absence of Ly a emission associated with the 7205 A Ly (3 emission peak. The Ly a 
emission would be seen at 8539 A but is completely absent (see Fig. 5 in White et al.) The optical depth in Ly (3 is smaller 
than Lya, so it is expected to have more light leak through the Ly/3 forest. However, the z ~ 6 Ly/3 trough also suffers 
from overlying z <~ 5 Ly a forest absorption, which reduces the incidence of Ly /3-only emission. And transparent Ly (3 
windows are always expected to have associated Ly a peaks. 

The combined optical depth 8 at 7205 A is T a (z~5) + Tp(z~6) = 0.97. The absence of any Lya emission at 8539 A 
implies a 3a lower limit r Q (8539) > 3.3, corresponding to 7)3 (7205) > 0.53. We conclude that T a (z~5) < 0.44. 

As White et al. point out, there are very few windows in the z = 5 Lya forest with such low optical depths. The 
probability that a randomly selected wavelength near 7200 A will fall in such a window is only about 1%. That is the 

7 And even among the coauthors of this paper. 

8 Here we adopt 2 = 6.42 for the quasar, as derived by Willott, McLure & Jarvis (2003) and confirmed by the CO observations of Bcrtoldi 
ct al. (2003). We also use a ratio Tp/r a = 0.16, correcting the 15% error pointed out by Songaila (2004). The effect of these changes on the 
computed optical depths is small. 
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fundamental problem: the strong peak in the Ly/3 trough requires the alignment of two unusually transparent windows, 
one in the Ly/3 forest and one in the Lya forest. However improbable it appears, this is the only explanation we see for 
the observed spectrum. 

Here we are using the approach criticized by Songaila (2004), simply scaling the measured r a and Tp values by the ratio 
of the transition strengths. This does in fact overestimate the T a /Tp ratio when there is a distribution of optical depths 
within a flux resolution element. However, the effect discussed by Songaila always reduces the Lya optical depths, so a 
simple scaling of Tp gives a conservative lower limit to the true optical depth. Suppose that the true Lya optical depth 
in some region has a distribution P(t), normalized such that / P{r)dr = 1. Then the observed mean optical depths in 
Ly a and Ly /3 are: 



T a = - In 

and 

rp = - In 



f 

Jo 

f 

Jo 



P(t) exp(— t)g?t 



P{t) cxp(— /r)dr 



(1) 
(2) 



where / = 0.16 is the nominal ratio between the Lya and Ly/3 optical depths. 

The exact value of the optical depth ratio depends on the distribution P(t), which has been explored using physical 
models of the gas density distribution and ionization balance (Fan ct al. 2002, Lidz et al. 2002, Cen & McDonald 2002, 
and references therein). For this paper we prefer to keep the equation in terms of the optical depth distribution P{t) 
rather rely on a more complicated physical model of the IGM to clarify the connections to directly observable quantities. 
A simple P(t) example illustrates the size of the effect. If P(t) is a Gaussian with mean f and width a, 

T a =lJ CT fcx[-.A/2V] 1 
\erfcx[(a 2 -f)A/2a] j 
where erfcx(x) is the scaled complementary error function, 

erfcx(a;) = -= / e**'* ' dt . (4) 



I IK 

When f > cr 2 , 

T « - f - y . (5) 

The equation for Tp is identical but with f and a replaced by ff and fa, respectively. 

The optical depth distribution effect explains the overall ratio of the mean transmissions observed in the Ly a and Ly [3 
Gunn-Peterson troughs for SDSS J1148+5251. Using the spectrum in Figure 1 (White et al. 2003), the mean optical 
depth in the Ly/3 trough (6.1 < z < 6.32) is Tp = 3.1 ± 0.1 while the mean in the Lya trough is r Q = 6.4 ± 0.3, which is 
smaller by a factor of three than might be expected from the optical depth ratio / = 0.16. (Note that in both absorption 
troughs in SDSS Jl 148+5251 we detect residual light, whereas in SDSS J1030+0524 we measure only an upper limit to 
the IGM transmission.) The transparency of the Lya trough is due to the optical depth distribution in the trough. If we 
apply the Gaussian model (cqn. 3) we find f = 26.4 and a = 10. Of course the Gaussian model used for this analysis is 
not particularly realistic, but the results will be similar for other P(t) models. 

In this case the t g / T a ratio is 0.48, which is (as demanded by the observations) indeed three times larger than / = 0.16. 
But the change is almost entirely the result of a reduction in the effective r Q compared with the true mean optical depth 
f. This effect is easy to understand: r Q begins to saturate when the high optical depth regions effectively transmit no 
light, and the mean transmission becomes dominated by light in low opacity regions. 

This is an important point: an optical depth distribution makes the IGM more transparent in the Lya absorption 
trough. As the Lya optical depth rises with increasing redshift, this effect tends to flatten the dT a /dz slope. However, 
we find that r Q increases more steeply with z starting at z ~ 6; correcting for the P(t) effect makes this change even 
more dramatic. The inclusion of this effect thus strengthens the evidence for a significant increase in the neutral hydrogen 
fraction beyond redshift 6. This is inconsistent with Songaila's (2004) conclusion that the inclusion of an optical depth 
distribution weakens the argument for an abrupt change in the ionization level at z > 6. 

Another way of stating this conclusion is that Tp/ f is actually a less biased estimator of the mean neutral hydrogen 
fraction (the physical quantity of interest) than is T a . The interpretation of the observed Ly [3 transmission is far less 
reliant on the details of P(r) than is the interpretation of the Lya transmission. Small changes in P(t) at t < 1 can 
have a large effect on r a , but Tp is sensitive only to P(t) changes integrated over the much larger range r < 1// w 6. The 
observed Tp optical depths increase very dramatically with redshift for the known z > 6 quasars (White et al. 2003; Fan 
et al. 2004). We interpret this as powerful direct evidence in favor of a rapid increase in the hydrogen neutral fraction at 
those redshifts. 

Returning to the question of the 7205 A emission in the Ly/3 trough of SDSS J1148+5251, it is clear that that the P(t) 
effect makes the absence of an associated Ly a peak even more difficult to understand because it increases the strength 
expected for the Lya peak. For example, the specific Gaussian P(t) described above (f/a — 2.6) is marginally consistent 
with the absence of Lya at 8539 A only if there is practically no z = 5 Lya forest absorption at all, t q (z=5) < 0.25. It 
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is possible to rule out P(r) distributions with stronger tails near r = 0. The properties of this system can consequently 
put limits on the permitted optical depth distribution P(t) at velocity resolutions finer than those directly observed. 

Perhaps in this quasar the z = 6.025 IGM is very transparent in Ly/3 yet opaque in Lya because the damping wings 
of the (roughly) neutral gas on either side are optically thick in Lya but not in Ly/3. If the presence or absence of 
damping wings can be shown to play a crucial role in interpreting the Ly [3 "leaks" in this system, then we may be able 
to strengthen the case for (or against) a substantial neutral hydrogen fraction in the IGM at z > 6. 

There is at least one satisfying aspect of this "two window" model for the 7205 A peak. A plausible way to create a 
transparent region in the Lya forest is for a star-forming galaxy close to the line of sight at z — 4.925 to ionize a bubble 
in the Lya forest. The halo of that same galaxy probably hosts the gas cloud that produces the strong C IV absorption 
doublet seen at 9200 A (z = 4.943, White et al. 2003). Note that our observations do not rule out a faint galaxy roughly 
coincident with the quasar (§3.1) even though the light in the 7205 A emission peak certainly comes from the quasar and 
not a galaxy. The angular scale at z — 4.9 is 6.3 kpc/arcsec (using the WMAP cosmology), so the galaxy need not lie too 
close to the line of sight to the quasar. A single unusual object along the line of sight can thus explain both the hole in 
the Ly a forest absorption and the strong high-ionization metal absorption. 

Finally, since the Lya and Ly/3 troughs of SDSS J1148+5251 do show a low level of detectable light (which is perhaps 
still slightly contaminated by the galaxy), this spectrum is worthy of additional detailed physical modeling as has been 
applied to SDSS J1030+0524. And given the apparently unusual characteristics of this line of sight, additional high 
rcdshift quasars sampling the IGM absorption in other directions would of course be of great value. 
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